Astrophysics and Space Science manuscript No. 

(will be inserted by the editor) 



Markus Bottcher 

Modeling the Emission Processes in Blazars 



o 
o 

(N 

3 ■ 

< 

m 



> 

m : 

00 1 

o : 

O ■ 

^ : 

Oh; 

6 ' 
a: 



13 



Received: date / Accepted: date 

Abstract Blazars are the most violent steady/recurrent 
sources of high-energy gamma-ray emission in the known 
Universe. They are prominent emitters of electromag- 
netic radiation throughout the entire electromagnetic spec- 
trum. The observable radiation most likely originates in 
a relativistic jet oriented at a small angle with respect 
to the line of sight. This review starts out with a general 
overview of the phenomenology of blazars, including re- 
sults from a recent multiwavelength observing campaign 
on 3C279. Subsequently, issues of modeling broadband 
spectra will be discussed. Spectral information alone is 
not sufficient to distinguish between competing models 
and to constrain essential parameters, in particular re- 
lated to the primary particle acceleration and radiation 
mechanisms in the jet. Short-term spectral variability 
information may help to break such model degenera- 
cies, which will require snap-shot spectral information 
on intraday time scales, which may soon be achievable 
for many blazars even in the gamma-ray regime with 
the upcoming GLAST mission and current advances in 
Atmospheric Cherenkov Telescope technology. In addi- 
tion to pure leptonic and hadronic models of gamma-ray 
emission from blazars, leptonic/hadronic hybrid models 
are reviewed, and the recently developed hadronic syn- 
chrotron mirror model for TeV 7-ray flares which are not 
accompanied by simultaneous X-ray flares ( "orphan TeV 
flares") is revisited. 
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1 Introduction 

Blazars (BL Lac objects and 7-ray loud flat spectrum ra- 
dio quasars [FSRQs]) are the most extreme class of active 
galaxies known. They have been observed at all wave- 
lengths, from radio through very-high energy (VHE) 7- 
rays. 46 blazars have been identified with high confi- 
dence as sources of > 100 MeV emission detected by 
the EGRET telescope on board the Compton Gamma- 
Ray Observatory |33U55| . and about one dozen blazars 
have now been detected at VHE 7-rays (> 350 GeV) by 
ground-based atmospheric Cerenkov telescopes (ACTs). 
Many of the EGRET-detected 7-ray blazars appear to 
emit the bulk of their bolometric luminosity at 7-ray en- 
ergies. Blazars exhibit variability at all wavelengths on 
various time scales. Radio interferometry often reveals 
one-sided kpc-scale jets with apparent superluminal mo- 
tion. 



1.1 Spectral classification of blazars 

The broadband continuum spectra of blazars are dom- 
inated by non-thermal emission and consist of two dis- 
tinct, broad components: A low-energy component from 
radio through UV or X-rays, and a high-energy com- 
ponent from X-rays to 7-rays. A sequence of blazar sub- 
classes, from FSRQs to low-frequency peaked BL Lac ob- 
jects (LBLs) to high-frequency peaked BL Lacs (HBLs) 
can be defined through the peak frequencies and relative 
vF v peak fluxes, which also seem to be correlated with 
the bolometric luminosity [21] • The sequence FSRQ — > 
LBL — > HBL is characterized by increasing vF v peak fre- 
quencies, a decreasing dominance of the 7-ray flux over 
the low-frequency emission, and a decreasing bolomet- 
ric luminosity (see Fig. 1). LBLs are intermediate be- 
tween the FSRQs and the HBLs. The peak of their low- 
frequency component is located at IR or optical wave- 
lengths, their high-frequency component peaks at several 
GeV, and the 7-ray output is of the order of or slightly 
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Fig. 1 SEDs of 3C 279 \34\ > BL Lacertae ]15l . and Mrk 501 \6'Tf . For each object, two simultaneous broadband spectra at 
two different epochs are shown. The curves show model fits, using a leptonic jet model. 



higher than the level of the low-frequency emission. How- 
ever, the existence and physical significance of this blazar 
sequence has recently been questioned and the apparent 
sequence attributed to selection effects due to the use of 
flux-limited samples (for a recent review see, e.g., |65p. 



1.2 Spectral variability of blazars 

Fig. n already illustrates that in particular the high- 
energy emission from blazars can easily vary by more 
than an order of magnitude between different EGRET 
observing epochs 89,60,61 . However, high-energy vari- 
ability has been observed on much shorter time scales, 
in some cases less than an hour |26| . 

BL Lac objects occasionally exhibit X-ray variability 
patterns which can be characterized as spectral hystere- 
sis in hardness-intensity diagrams (e.g., [Hal3QI25l91| ). 
This has been interpreted as the synchrotron signature 
of gradual injection and/or acceleration of ultrarelativis- 
tic electrons in the emitting region, and subsequent ra- 
diative cooling (e.g., [39112 7U36U43U45ll9| ) . While spectral 
hysteresis has so far only been clearly identified in HBLs, 
it should also occur in the soft X-ray emission of LBLs 
if their synchrotron component extends into the soft X- 
ray regime. However, LBLs are generally fainter at X- 
ray energies than HBLs, making the extraction of time- 
resolved spectral information observationally very chal- 
lenging. Fig. [2i shows the results of a BeppoSAX obser- 
vation of BL Lacertae in 2000 |71U15| . Rapid flux and 
spectral variability of blazars is also commonly observed 
in the optical regime, often characterized by a spectral 
hardening during flares (see, e.g., Fig. or |4*5I57) 1. 



2 Preliminary results from 3C279 2006 

The flat spectrum radio quasar 3C279 (z — 0.538) is one 
of the best observed flat spectrum radio quasars, not 
at last because of its prominent 7-ray flare shortly af- 
ter the launch of the Compton Gamma-Ray Observatory 
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Fig. 3 Preliminary optical (BVRI) light curves of 3C 279 
over the entire core period of the 2006 multiwavelength cam- 
paign. The shaded area indicates the period of the INTE- 
GRAL X-ray and soft 7-ray observations. The red horizontal 
line indicates the trigger criterion of R = 14.5. 



(CGRO) in 1991. It has been persistently detected by the 
Energetic Gamma-Ray Experiment Telescope (EGRET) 
on board CGRO each time it was observed, even in its 
very low quiescent states, e.g., in the winter of 1992 - 
1993, and is known to vary in 7-ray flux by roughly two 
orders of magnitude 51,90 . It has been monitored in- 
tensively at radio, optical, and more recently also X-ray 
frequencies, and has been the subject of intensive mul- 
tiwavelength campaigns (e.g., [HOI)- Its broadband SEDs 
at several epochs are rather well determined, but a com- 
plete compilation and modeling (using a leptonic model, 
see t|3.1fl of all available SEDs simultaneous with the 11 
EGRET observing epochs yielded somewhat inconclu- 
sive results |34) . Furthermore, in spite of the intensive 
past observational efforts, the broadband spectral vari- 
ability of 3C 279 is still rather poorly understood (see, 
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Fig. 2 Spectral variability fitting of BL Lacertae in 2000 |1_JJ: (a) Time-dependent leptonic fits to the Nov. f high state; 
(b) Various hadronic fits, differing mainly in their co-moving magnetic-field and synchrotron photon energy density, (c) 
Comparison of the simulation results corresponding to fits from panel (a) and (b) to the optical color-magnitude correlation; 
(d) comparison of the leptonic fit results to the X-ray hardness-intensity diagrams during a short flare observed by BeppoSAX 
on Nov. f . Based on count rates in the three BeppoSAX NFI energy channels LECS [0.5 - 2 keV], MECS [2-4 keV], and 
MECS [4 - fO keV], and X-ray hardness ratios: HR1 = MECS [2-4] /LECS [0.5-2], HR2 = MECS [4- fO] /MECS [2-4]. 



For the reasons stated above, we |20] proposed an 
intensive multiwavelength campaign in an optical high 
state of the 3C 279, in order to investigate its corre- 
lated radio - IR - optical - X-ray - soft 7-ray variabil- 
ity. The campaign was triggered on Jan. 5, 2006, when 
the source exceeded an R-band flux corresponding to R 
= 14.5. It involved intensive radio, near-IR (JHK), and 
optical monitoring by the WEBT collaboration through 
March of 2006, focusing on a core period of Jan. and 
Feb. 2006. X-ray and soft 7 observations were carried out 
by all instruments on board the International Gamma- 
Ray Astrophysics Laboratory (INTEGRAL) during the 
period of Jan. 13 - 20, 2006. Additional, simultaneous 
X-ray coverage was obtained by Chandra and Swift XRT. 
These observations were supplemented by extended X- 
ray monitoring with the Rossi X-Ray Timing Explorer 



(RXTE) and VLBA monitoring at 43 GHz. The analysis 
of the data collected during this campaign is currently in 
progress. Here, first, preliminary results from this cam- 
paign are presented. Final analysis results of the WEBT 
(radio - IR - optical) campaign will be published in [2j , 
while a comprehensive report on the result of the entire 
multiwavelength campaign will appear in |20j . 

Fig. shows the preliminary optical (BVRI) light 
curves, including about 3/4 of all collected data from the 
campaign. The figure illustrates that the source showed 
substantial, closely correlated variability in all wave bands 
throughout the entire core campaign. Between Jan. 8 and 
15 (i.e., including the time of the INTEGRAL observa- 
tions), the optical flux was persistently fading in all opti- 
cal bands so that the R band flux was actually below the 
intended trigger threshold. However, this may not nec- 



4 



Markus Bottcher 



3C279 Multiwavelength Lightcurves 



Swift XRT (0.1 cts/s) 



INTEGRAL IBIS(10 cm s ) 



5 GHz 
8 GHz 

22 GHz 
37 GHz 



i f 



746 748 750 752 754 
JD - 2,453,000 



756 



3C279 



Fig. 4 Preliminary X-ray and soft 7-ray (top panel), opti- 
cal (middle panel), and radio (bottom panel) light curves of 
3C 279 during the time of the INTEGRAL observations on 
Jan. 13 - 20, 2006. 
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Fig. 5 Simultaneous spectral energy distributions of 3C 279 
during 4 observing epochs, including a snapshot SED on Jan. 
15, 2006 (red triangles pointing down). 



essarily be bad news since this might allow us to study 
the correlated radio through X-ray flux decay time scales 
during a period of clean, steady decay after a major op- 
tical outburst. We might thus be able to probe energy- 
dependent electron cooling time scales throughout the 
entire synchrotron component and the low-energy part 
of the high-energy component (generally attributed to 
synchrotron-self-Compton emission in leptonic jet mod- 
els, e.g., P2) of the SED of 3C279. Fig.Hillustrates that 
the general trend of flux decline throughout most of the 
time window of the INTEGRAL observations was also 
shared by the high-frequency (37 GHz) radio flux, while 
the X-ray flux detected by Swift indicates a slowly ris- 
ing trend after Jan. 14. Possible interpretations of this 
behavior will be discussed after the brief introduction to 
leptonic blazar jet models in W3.l\ 

In Fig. [31 we compare a snapshot SED of 3C 279 
of Jan. 15, 2006 to several other simultaneous SEDs, in- 
cluding the exceptional June 1991 flare and the very low, 
quiescent state of Dec. 1992 - Jan. 1993, as well as the 
time of the previous INTEGRAL observation in June 
2003 [TJJ. It reveals a surprising result: Even though the 
optical flux is well (factor 2-3) above the typical low- 
state values, the X-ray and soft 7-ray fluxes and spec- 
tra are in the range of the lowest X-ray fluxes measured 
from the source, including the 1992/93 and 2003 quies- 
cent states. Before we procede with a possible physical 
interpretation of the results from this campaign, let us 
briefly review the salient features of blazar jet models. 



3 Blazar models 

The high inferred bolometric luminosities, rapid vari- 
ability, and apparent superluminal motions provide com- 
pelling evidence that the nonthermal continuum emission 



of blazars is produced in <, 1 light day sized emission 
regions, propagating relativistically along a jet directed 
at a small angle with respect to our line of sight. It is 
generally agreed that the low-frequency component of 
blazar SEDs might be synchrotron radiation from non- 
thermal, ultrarelativistic electrons. Several electron in- 
jection/acceleration scenarios have been proposed, e.g. 
impulsive injection near the base of the jet (e.g., [221E3I; 
such a scenario might also apply to originally Poynting- 
flux dominated jets, see (EH)) isolated shocks propagat- 
ing along the jet (e.g., |S^ll39ll78ll5T] '). internal shocks 
from the collisions of multiple shells of material in the 
jet [531, stochastic particle acceleration in shear bound- 
ary layers of relativistic jets (e.g., H2][73]), magnetic re- 
connection in Poynting-flux dominated jets (e.g., |80| ) . 
or hadronically initiated pair avalanches |37| . 

Significant progress has recently been made in our un- 
derstanding of particle acceleration at relativistic shocks 
(e.g., 1,64 ,62 88 ) and the conversion of relativistic bulk 
kinetic energy into relativistic particles and ultimately 
into radiation 69,76,86,53 . While particle acceleration 
at relativistic parallel shocks might produce electron in- 
jection spectra of Neil) oc j~ q with 2.2<q<2.3 (e.g., 
PP), oblique shocks tend to produce much softer injec- 
tion spectral indices (e.g., EU). In contrast, the effect 
of stochastic acceleration in resonance with plasma wave 
turbulence behind the shock front might harden the in- 
jection index significantly, possibly even beyond q ~ 1 
|88| . In the framework of an inhomogeneous jet with 
a fast inner spine and a slower, outer cocoon, particle 
acceleration at shear boundary layers may become the 
dominant acceleration mechanism | 63U84U73| . This may 
lead to the development of relativistic particle spectra 
with indices of q < 2 and a high-energy bump around 
the energy where the particle acceleration rate equals 
the energy loss rate. 
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However, the lack of knowledge of the primary jet 
launching mechanism (poynting-flux dominated vs. mag- 
neto- hydrodynamic, see, e.g., [SOI ) an d the primary pair 
loading of the jet are currently severe problems in ty- 
ing the properties of the particle acceleration mecha- 
nism to more fundamental physical properties of the 
accretion-powered disk-jet system. There is mounting ev- 
idence that — if the high-energy emission of blazars is 
leptonically dominated — jets of blazars might be en- 
ergetically and dynamically dominated by their proton 
content, though pairs may still greatly outnumber pro- 
tons 



While the electron-synchrotron origin of the low-fre- 
quency emission is well established, there are two funda- 
mentally different approaches concerning the high-energy 
emission. If protons are not accelerated to sufficiently 
high energies to reach the threshold for p^j pion produc- 
tion on synchrotron and/or external photons, the high- 
energy radiation will be dominated by emission from ul- 
trarelativistic electrons and/or pairs (leptonic models). 
In the opposite case, the high-energy emission will be 
dominated by cascades initiated by p~f pair and pion 
production as well as proton, % , and ^ synchrotron 
radiation (hadronic models). 

3.1 Leptonic blazar models 

If protons are not accelerated to sufficiently high ener- 
gies to reach the threshold for pr/ pion production on 
synchrotron and/or external photons, the high-energy 
emission will be dominated by ultrarelativistic electrons 
and/or pairs. In this case, high-energy emission can be 
produced via Compton scattering off the same ultrarela- 
tivistic electrons producing the synchrotron emission at 
lower frequencies. Possible target photon fields are the 
synchrotron photons produced within the jet (the SSC 
process: |52U5()U8p or external photons (the EC process). 
Possible sources of external seed photons are accretion- 
disk photons entering the jet directly or after be- 
ing reprocessed in the circumnuclear material (e.g., the 
broad line regions of quasars) |ZZK3, jet synchrotron 
emission reflected off clouds in the circumnuclear mate- 
rial |31| . infrared emission from a dust torus around the 
central engine 6 3 . In addition, 77 absorption, pair pro- 
duction, and synchrotron self absorption must be taken 
into account in a self-consistent leptonic blazar model. As 
the emission region is propagating relativistically along 
the jet, continuous particle injection and/or acceleration 
and subsequent radiative and adiabatic cooling, particle 
escape, and possibly also the deceleration of the jet, in 
particular in HBLs 28 32 , have to be considered. 

As one may expect from the basic description in the 
previous paragraph, leptonic models (especially when con- 
sidering a substantial contribution from external radia- 
tion fields as targets for Compton upscattering to 7-rays) 
require the specification of a rather large number of pa- 
rameters. Several of these parameters can be estimated 



to a reasonable degree of accuracy from observables (see, 
e.g., ^H] for a discussion of such parameter estimates in 
the case of BL Lacertae during the multiwavclength cam- 
paign of 2000). In particular, there are several ways to 
obtain an estimate on the magnetic field in the emission 
region, which is an essential ingredient to also get a han- 
dle on the efficiency of Fermi and stochastic acceleration 
of particles in the jet. 

One way to estimate the co-moving magnetic field 
can be found by assuming that the dominant portion of 
the time-averaged synchrotron spectrum is emitted by a 
quasi-equilibrium power-law spectrum of electrons with 
-^3(7) = n o Vb 7 _p for 71 < 7 < 72; here, Vb is the 
co-moving blob volume. The normalization constant no 
is related to the magnetic field through an equiparti- 
tion parameter cb = UB/u e (in the co- moving frame). 
Note that this equipartition parameter only refers to 
the energy density of the electrons, not accounting for a 
(possibly greatly dominant) energy content of a hadronic 
matter component in the jet. Under these assumptions, 
the vF v peak synchrotron flux /| y at the dimensionless 
synchrotron peak energy e sy = his sy / (m e c 2 ) is approxi- 
mately given by 

f? = ( D B ) V2 £££([1 + *] ^ ^ CD 



288 d\ 



where D = 10 L>i = (r[l - j3 r cos o b s ]) _1 is the Doppler 
boosting factor, = 10 27 d%j cm is the the luminos- 
ity distance of the source, and B CI = 4.414 x 10 13 G. 
The electron spectrum normalization used to derive Eq. 
^ is based on the presence of a power-law shape with 
a photon energy index a > 1 (with F v oc u~ a ) of the 
synchrotron spectrum beyond the synchrotron peak. If 
the SED reveals such a spectral shape, the underlying 
electron spectrum always has an index of p > 3. Eq. 
then yields a magnetic-field estimate of 



B, 



9£>r 



d 4 f 2 p 2 
"27 J-10 e B 



[l + z] 4 e sy ^ 6 R%[p-2] 



1/7 



G. 



(2) 



where /_ 10 = f?/(lQ- ia ergs cm" 2 s^ 1 ), e sy ,_ 6 = e sy /l0' 
and Rb = 10 15 -R15 cm is the transverse radius of the 
emission region. The bulk Lorentz factor r and Doppler 
factor D can usually be constrained from superluminal 
motion measurements and/or from constraints on the 
compactness of the emission region, which is an approxi- 
mate measure of the optical depth to 77 absorption. An 
estimate of the size of the emission region can be inferred 
from the minimum variability time scale, i var = t vai! h hr, 
as Rb ~ 10 15 D\ i V ar,h cm. 

We can apply this estimate to the results of our 2006 
campaign on 3C 279 described in £f2I For z = 0.538, 
H = 70 km s- 1 Mpc" 1 , Q m = 0.3, and Q A = 0.7, the 
luminosity distance is dh = 9.3 x 10 27 cm. From previous 
work on 3C 279 (see, e.g., for a summary) one finds 
typical values of D ~ 10, and Rb ~6x 10 16 cm. Unfor- 
tunately, the synchrotron peak seems to lie in the mm 



6 



Markus Bottcher 



- far-IR regime of the spectrum, which was not covered 
during the campaign. Thus, its position is not very well 
constrained. Visual inspection of the SED suggest values 
of /-io ~ 0.7 and v sy w 8 x 10 13 Hz, corresponding to 
e sy ,_ 6 w 0.65. The IR - optical (UBVRIJHK) spectral 
index is a — 1.75, which corresponds to p — 4.5. These 
values lead to a magnetic-field estimate of 

B eB ^0.6D^e 2 J 7 G. (3) 

which appears to be relatively weak compared to typi- 
cal values of ~ a few G inferred from modeling efforts 
on other FSRQs. The steep spectral index suggests that 
electrons emitting IR and optical radiation might al- 
ready be in the strong cooling regime, where the ra- 
diative cooling time is shorter than the dynamical time 
scale. Even in that case, one would infer that electrons 
are injected through the primary particle acceleration 
mechanism with an injection index of q ~ 3.5, which 
would imply rather inefficient acceleration, and may point 
towards Fermi acceleration at oblique shocks as the pri- 
mary particle acceleration mechanism (see the discussion 
in 

Although several model parameters of blazar jet mod- 
els can be reasonably well constrained from the broad- 
band spectral properties of blazars, spectral fitting alone 
is generally insufficient to constrain all relevant model 
parameters (see, e.g., the analysis in ^H] for the case of 
the LBL W Comae). Thus, it is now widely agreed that 
spectral and variability properties of blazars have to be 
taken into account simultaneously in order to extract as 
much physical information as possible from simultaneous 
multiwavelength observing campaigns. 

Significant progress has been made in the past few 
years to combine spectral and variability modeling of 
blazars using leptonic models. In particular, the spectral 
variability of HBLs has been modelled in great detail 
by many authors, using pure SSC models (e.g., [231211 
I4UU81| . Time-dependent blazar modeling including ex- 
ternal soft photon sources is also advancing rapidly (e.g., 
|78ll9"ll44U82| 'l. providing tools for the interpretation of si- 
multaneous spectral and variability data from LBLs and 
FSRQs. An instructive example of combined fitting of 
SEDs and rapid spectral variability of BL Lacertae can 
be found in (see also Fig. EJ. 

In this context, one should mention an alternative 
way of estimating the magnetic field in blazar jets. This 
is based on a possible time delay between light curves at 
two different frequencies at which the emission is dom- 
inated by synchrotron emission. Assuming that such a 
delay is caused by synchrotron cooling of high-energy 
electrons with characteristic observed synchrotron pho- 
ton energy -E S y,o = Eq keV to lower energies with corre- 
sponding synchrotron energy E syi x = E\ keV, the mag- 
netic field can be estimated as: 

B dclay = 0AD; 1/3 (1 + fc)- 2 / 3 (At^T 2/3 



(E; 1/2 - E- 1/2 f/3 G. (4) 

where k = Mph/^B is the ratio of energy densities in the 
photon field in the frame co-moving with the emission 
region and the magnetic field, and Z\i£ bs is the observed 
time delay in hours. 

Our preliminary analysis of the results from the 2006 
campaign on 3C 279 did not yield any evidence for time 
lags between any of the near-IR and optical bands. Fur- 
thermore, the SED of Jan. 15 (see Fig. EJ) clearly shows 
that the X-ray emission is dominated by the low-energy 
end of the high-energy spectral component. Therefore, 
possible X-ray — optical delays can not be used for the 
estimate in Eq. 0J 

The somewhat surprising apparent trend of a slow 
rise in the X-ray and soft 7-ray fluxes during the con- 
tinuing optical fading trend illustrated in Fig. 0] may be 
interpreted in the following way: In leptonic jet models, 
the low-energy end of the high-frequency bump of the 
SEDs of blazars (which covers the X-ray band in the case 
of 3C 279) is generally attributed to synchrotron self- 
Compton emission from rather low-energy electrons. The 
differential number density of electrons at moderately 
relativistic energies will be gradually built up through 
the relatively long electron cooling time scale of elec- 
trons at those energies. In addition, the self-generated 
synchrotron photon field, serving as a target for the SSC 
process will also be gradually built up throughout at 
least the dynamical time scale. Consequently, substantial 
time delays of the X-ray emission behind optical flaring 
activity may result (see, e.g., |7S] for a detailed discus- 
sion). This could mean that the slow rise of the X-ray 
flux is, in fact, the delayed response to the optical flare 
observed around Jan. 8, 2006. A more quantitative anal- 
ysis of this interpretation will be presented in |2"U] . 

3.2 Hadronic blazar models 

If a significant fraction of the kinetic power in the jet is 
converted into the acceleration of relativistic protons and 
those protons reach the threshold for pj pion produc- 
tion, synchrotron-supported pair cascades will develop 
[491148] . The acceleration of protons to the necessary ul- 
trarelativistic energies requires high magnetic fields of at 
least several tens of Gauss. In the presence of such high 
magnetic fields, the synchrotron radiation of the primary 
protons [2"ll57| and of secondary muons and mesons |70l 
l5Tll58ll59j must be taken into account in order to con- 
struct a self-consistent synchrotron-proton blazar (SPB) 
model. Electromagnetic cascades can be initiated by pho- 
tons from 7r°-decay ("7r° cascade"), electrons from the 
7T^~ — > /i ± — » e ± decay ("7r ± cascade"), p-synchrotron 
photons ( "p-synchrotron cascade"), and fi-, n- and K- 
synchrotron photons ("/^-synchrotron cascade"). 

It has been shown in [581159] that the "ir° cascades" 
and "7^ cascades" generate featureless 7-ray spectra, 
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in contrast to "p-synchrotron cascades" and "/i -syn- 
chrotron cascades" that produce a double-bumped 7-ray 
spectrum. In general, direct proton and ^ synchrotron 
radiation is mainly responsible for the high energy bump 
in blazars, whereas the low energy bump is dominanted 
by synchrotron radiation from the primary e~, with a 
contribution from secondary electrons. Fig.JSJ} shows fits 
to the SED of BL Lacertae in 2000, using the hadronic 
SPB model |55] . 

3.3 Blazar unification 

Leptonic models have been used successfully to repro- 
duce simultaneous SEDs of several blazars. Spectral mod- 
eling results are now converging towards a rather consis- 
tent picture The sequence HBL -> LBL -> FSRQ 
appears to be related to an increasing external-Compton 
contribution to the 7-ray spectrum. While most FSRQs 
are successfully modelled with EC models (e.g., [23TT7H 
E2EH), the SEDs of HBLs are consistent with pure SSC 
models (e.g., 54, 68 , 67 40 ) . LBLs (e.g., BL Lacertae, 
W Comae) often seem to require an EC component to 
explain their EGRET spectra f751i47lllllllO| . One gen- 
erally finds that HBLs require higher average electron 
energies and lower magnetic fields than LBLs and FS- 
RQs. In the framework of a unified leptonic model, this 
basic parameter sequence may be related to an increas- 
ing importance of EC cooling along the sequence HBL 
-> LBL -> FSRQ [2S]- It has been suggested that the de- 
creasing importance of external radiation fields along the 
sequence FSRQ — > LBL — » HBL may be an evolution- 
ary effect related to the gradual depletion of a limited 
reservoir of circumnuclear material [2111181113) . 

Hadronic blazar models also offer a physical inter- 
pretation for the spectral sequence of BL Lac subclasses 
|59| . The spectra of HBLs are well reproduced by p- 
synchrotron dominated SPB models where the intrin- 
sic primary synchrotron photon energy density is small, 
consistent with the low bolometric luminosity of those 
objects. As the synchrotron photon energy density in- 
creases towards LBL-like synchrotron properties, pro- 
tons suffer increasingly strong pj pion production losses, 
and the contributions from the and ^ synchrotron 
cascades become increasingly dominant at higher ener- 
gies. This results in a decreasing vF v peak frequency 
of the 7-ray component. The effect of external photon 
sources might further enhance the ^-synchrotron and 
/^-synchrotron cascade contributions, reproducing the 
transition to quasar-like properties. 

It should be pointed out here that the blazar se- 
quence, if real, can be explained in the framework of 
both leptonic and hadronic models, but it is not a predic- 
tion of either one of the classes of models. Consequently, 
even if future observations reveal evidence for a rather 
uniform distribution of peak frequencies and relative lu- 
minosities between the two main spectral components of 
blazars, both model variants remain generally viable. 



3.4 Hybrid blazar models 

The leptonic and hadronic models discussed above are 
certainly only to be regarded as extreme idealizations of 
a blazar jet. Realistically, both types of processes might 
play a role to some extent and should thus be considered 
to a comparable level of sophistication. The recent obser- 
vation of isolated TeV flares without simultaneous X-ray 
flares (a phenomenon sometimes referred to as "orphan 
TeV flares") in 1ES 1959+650 [HI and Mrk 421 \7\ may, 
in fact, provide rather strong support for the importance 
of hadronic processes in objects of which other spectral 
and variability features arc generally well reproduced by 
leptonic jet models, since standard leptonic SSC models 
predict a close temporal flux correlation between the syn- 
chrotron and Compton components. In 1ES 1959+650, 
the "orphan" TeV flare was preceded by an ordinary, 
correlated X-ray and TeV-flare, which can be generally 
well understood in the context of leptonic SSC models. 
This finding strongly suggests the need for models that 
explain flares dominated by leptonic interactions as well 
as flares where non-leptonic components might play an 
important role within the same system. 

Hadronic processes in the context of models with lep- 
tonically dominated blazar emission have been consid- 
ered by several authors, e.g.: 

— A "supercritical pile" model was suggested in |37j. 
In this model a runaway pair production avalanche 
is initiated by mildly relativistic protons interact- 
ing with reflected synchrotron photons via pry pair 
production, as the primary pair injection mechanism 
in blazar jets. Spectral characteristics resulting from 
this model as applied to gamma-ray bursts have been 
considered in [5"3"|. 

— The conversion of ultrarelativistic protons into neu- 
trons via 737 pion production on external soft photons 
was suggested in 0] as a mechanism to overcome syn- 
chrotron losses of protons near the base of blazar jets 
and, thus, to allow blazar jets to remain collimated 
out to kpc scales. 

— Focusing on applications to Gamma-ray bursts, in 
|66j a fully self-consistent, time-dependent homoge- 
neous one-zone model was developed for the radiation 
from a relativistic plasma which assumes electron and 
proton injection into a power-law distribution and 
includes the self-consistent cooling of protons by p-f 
pion production processes and their contributions to 
the pair populations (and their radiative output). 

— In the hadronic synchrotron mirror model |16U72| , de- 
veloped specifically to explain the "orphan" TeV flare 
phenomenon in 1ES 1959+650 mentioned above, the 
primary, correlated X-ray and TeV flare is explained 
by a standard SSC model while the secondary TeV- 
flare is explained by 7r°-decay 7-rays as a result of 
photomeson production from relativistic protons in- 
teracting with synchrotron photons that have been 
reflected off clouds located at pc-scale distances from 
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the central engine. Since some of the numerical values 
in the original paper ^H] were in error (see the erra- 
tum) , this model will be revisited and re-evaluated in 
the next section. 



4 The hadronic synchrotron mirror model 
revisited 

The recent "orphan" TeV flare of 1ES 1959+650 led 
to the development of the hadronic synchrotron mirror 
model [ToTr72| . The basic model geometry is sketched in 
Fig. H3 A blob filled with ultrarelativistic electrons and 
relativistic protons is traveling along the relativistic jet, 
defining the positive z axis. Particles are primarily accel- 
erated very close to the central engine (Fl) in an explo- 
sive event which is producing the initial synchrotron + 
TeV flare via the leptonic SSC mechanism. Acceleration 
of relativistic particles is expected to persist throughout 
the further propagation of the emission region along the 
jet through any of the processes (internal shocks, shear 
layer boundary acceleration, etc.) mentioned in thus 
leading to a sustained level of (quiescent state) UV/X- 
ray synchrotron emission, as observed in 1ES 1959+650 
|4T] . A fraction of this synchrotron radiation is reflected 
off a gas cloud (the mirror) located at a distance R m from 
the central engine. For the sake of analytical tractability, 
we assume that the mirror (M) is a homogeneous shell 
with a reprocessing depth r rn = 0.1 r_i. 

The observed time delay between the primary syn- 
chrotron flare and the secondary flare due to interactions 
of the blob with the first reflected synchrotron flare pho- 
tons to re-enter the blob was Z\£ bs = 20 At2o days, and is 
related to the distance of the reflector by At b s ~ 2T^c - 
Thus, R m w 3r , 1 2 Z\t 2 o pc As shown in [TJ], a cloud of 
reflecting gas of typical radial extent Ar = 10 17 Ar\i cm 
and density n c = 10 6 ng cm~ 3 at this distance from the 
central source will result in a negligible flux in emission 
lines, and the expected duration of the secondary flare 
is wg bs ~ 1.2 Pj~ hr, consistent with the observed time 
profile of the secondary TeV flare in 1ES 1959+650. 

From the observed vF v fluxes in synchrotron and 
TeV emission during the secondary TeV flare, vF v (sy) ~ 
2 x 10" 10 ergs s" 1 cm" 2 and vF v (600 GeV) - 3 x 10" 10 
ergs s _1 cm~ 2 [41 j . we find the co-moving luminosities, 
L' sy ~ 1.0 x lO 41 ^ -4 ergs s" 1 and L' VHE ~ 1.5 x 
10 41 Tf 4 ergs s -1 . With a blob radius of R' B = 10 16 R 16 
cm, this yields a co-moving synchrotron radiation field 
of 



9dl 



sy 



AR'%cr A 



vF v {sy) 



6.0 x 10 



-3 



.Tf 4 i?^ 2 ergs cm" 3 . 



(5) 



The characteristics of the reflected synchrotron flux (with 
co- moving energy density u'ruy) have been calculated in 



, Secondary Flare 




pos. synchrotron 
(opt. - UV - soft X-ray) 



Primary Flare: 
el. synchrotron + SSC 

Fig. 6 Geometry of the model. A primary synchrotron flare 
is produced by the emission region near the center of the sys- 
tem (Fl). Synchrotron emission is reflected at the mirror (M), 
and re-enters the emission region. A quiescent level of syn- 
chrotron emission towards the mirror will be sustained follow- 
ing the initial synchrotron flare. Its reflection into the emis- 
sion region is the primary source of external photons leading 
to the secondary, "orphan" TeV flare as the emission region 
approaches the mirror (F2). 



[T2| . Using their Eq. (4) in the limit z > 2rR' B , we 
find 



u' Rsy ~ 0.24 1^ 1 T- 1 (Ar^ 1 ) ergs cm 



(6) 



The co-moving luminosity from pj — > A — > p + ir° — * 
p + 2j produced by protons of a given energy 7' is then 
given by 

ti 8 1 1 a 1 70MeV . , . 

i'vHE -ca 4 uW/p Ai p — N' p {i p ). (7) 

<J Rsy 

where a a ~ 300 jA> is the A resonance peak cross sec- 
tion, Ay' p ~ Yp/2 parametrizes the FWHM of the A 
resonance, and N' ' p (pi' „) is the differential number of pro- 
tons at energy 7' . With this, the observable vF v peak 
flux in the TeV flare can be estimated as 



vF u (VBE) 



VHE 



~3.6x 10- 59 V' p (7' P )i? sy 3 1 

r_i (Am)- 1 ergs cm" 2 s" 1 . (8) 

Setting this equal to the observed VHE peak flux, we 
find 



N' p (i p ) ~ 8.3 x 10 48 rZ\ Rie Ari7. 



(9) 



Modeling the Emission Processes in Blazars 
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If the non-thermal proton spectrum is a straight power- 
law with energy index s = 2 and low-energy cut-off 
t'p min = r> this corresponds to a total proton num- 
ber of N' p . total ~ 7.5 x 10 55 rf 3 rli 1 £ sy ,i #i 6 and 
a proton number density of 

n' p ~ 1.8 x 10 7 Tf 3 E sy>1 r^ 1 i?^ 2 ^n 7 cm" 3 . (10) 

We note that, in order to bring this to a value in the 
range of electron densities typically invoked for leptonic 
jet models of blazar emission (n' e ~ a few 10 3 - 10 4 ), a 
substantially higher Doppler factor, a much more com- 
pact mirror, or a flatter proton spectrum, s < 2, seems 
to be required. From Eq. we find the total energy 
in relativistic protons in the blob as E'b. p ~ 8.5 x 
1Q 47 T^ 2 E syA Til R 16 Ar 17 er S- If blobs of such rela " 
tivistic plasma fill a fraction / = 10~ 3 /_3 of the jet, this 
corresponds to a kinetic luminosity in protons of L p ~ 
1.8 x 10 48 /_3 E sy ^i tZ\ Ar\7 ergs s _1 . These numbers in- 
dicate that orphan TeV flares as observed in 1ES 1959+650 
and Mrk 410 seem to require some rather extreme condi- 
tions which might only be present at very rare occasions. 
This may explain why there have so far only been very 
few examples observed. 
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